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(■ — ' Abstract. A programme is outlined for the assembly of a comprehensive dielectronic recombination database within the gen- 

, eralized collisional-radiative (GCR) framework. It is valid for modelling ions of elements in dynamic finite-density plasmas 

■ such as occur in transient astrophysical plasmas such as solar flares and in the divertors and high transport regions of magnetic 

fusion devices. The resolution and precision of the data are tuned to spectral analysis and so are sufficient for prediction of the 
dielectronic recombination contributions to individual spectral line emissivities. The fundamental data are structured according 
to the format prescriptions of the Atomic Data and Analysis Structure (ADAS) and the production of relevant GCR derived data 
for application is described and implemented following ADAS. The requirements on the dielectronic recombination database 
are reviewed and the new data are placed in context and evaluated with respect to older and more approximate treatments. 
Illustrative results validate the new high-resolution zero-density dielectronic recombination data in comparison with measure- 
ments made in heavy-ion storage rings utilizing an electron cooler. We also exemplify the role of the dielectronic data on GCR 
coefficient behaviour for some representative light and medium weight elements. 
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5_j 1. Introduction ate autoionizing states - indeed, entire Rydberg series. It was 

the significance of the effective statistical weight of so many 
Dielectronic recombination (DR) is the dominant electron-ion states ^ led Burgess (1%4) to recogn i Z e the im- 

recombination process in many astrophysical and laboratory portance of dielectronic rec0 mbination in the first instance. A 

plasmas. It plays an important role in determining both the consequence of this is mat a huge popu i ation structure mode l 

level populations and the ionization balance of both high- and ifj required fa prindple foj . further progress In astrophysics 

low-temperature non-LTE plasmas over a wide range of elec- and ^ problem hag usuaUy been mad£ manageable 

tron densities. Dielectronic recombination can be viewed as by simply summing over all final states so as t0 produce a t0 . 

a two-step process. Firstly, a free electron excites an electron ^ dielectronic reC ombination rate coefficient. In combination 

in an ion X\ say, and in the process transfers sufficient of wkh mdiati ve recombination, this is the effective recombina- 

its energy that it is captured into an autoionizing state of the tion coefficient from me point of view of an ionization state 

ion X*-\ The process is reversible since the total energy of Qnly _ ±e so . called l(;oronal picture - In the latter> e xcited-state 

the system remains conserved. However, if an electron (either populations are depleted exclusively by spontaneous radiative 
the captured electron or an electron in the parent core) then and are small compared to tho se of ground states 

makes a spontaneous radiative transition that leaves the ion in _ wkh which they ^ ^ quasi . static equilibrium. Collisional 

a non-autoionizing (bound) state then the recombination can be processes are negligible, except with ground state targets. In 

viewed as complete, at least in the low-density (coronal) limit. ^ ±Q ionization state is determined by balancing the effec- 

Dielectronic recombination can take place via many intermedi- tiye recombination rate of (the ground state of) x+ , against the 

~ , ~ . ~„ „ . „ collisional ionization rate of (the ground state of) X +z_1 . 
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For dynamic finite-density plasmas, there are therefore two 
critical limitations to the coronal approximation which must 
be addressed. Both require a re-appraisal of the strategy for 
computing dielectronic recombination, which is the objective 
of this paper. Firstly, from the atomic point of view, a dynamic 
plasma is one for which the timescale of change in plasma 
parameters (especially electron temperature, T e , and electron 
density, N e ) is comparable with the lifetime of metastable pop- 
ulations of its constituent ions. This situation is a concern, for 
example, for impurities near plasma contacted surfaces in fu- 
sion (Summers et al. 2002) and for active solar events (Lanza et 
al. 2001). The population of an ionization stage may no longer 
be assumed to be concentrated in the ground state of each ion. 
A significant population is found in the metastables, and they 
are not in quasi-static equilibrium with the ground state. This 
means that such metastables may be the starting point for re- 
combination events, and so the time-evolution of their popu- 
lations must be tracked in the same manner as for the ground 
states. We call this the generalized collisional-radiative (GCR) 
picture (Summers & Hooper 1983). 

The second issue is that of finite-density effects. The coro- 
nal (zero-density limit) picture assumes that, following the two- 
step dielectronic recombination process, the resultant (non- 
autoionizing) excited-state (electron) radiatively cascades back 
down to the ground state without collisional disruption. At fi- 
nite electron densities, this radiative cascade can be interrupted 
by further electron collisions which redistribute the popula- 
tion - in particular, ionization (possibly stepwise) out of ex- 
cited states, which reduces the effective dielectronic recombi- 
nation rate. Collisional-radiative modelling removes the lim- 
itations of the coronal model, but at the cost of much more 
elaborate excited population calculations. These in turn require 
much more detailed dielectronic recombination data, and in an 
easily accessible form. The onset of these density effects on 
the 'post-DR' population structure depends markedly on ion 
charge, but can be significant even at electron densities as low 
as N e Z 10 8 ctrT 3 typical of the solar corona. At much higher 
densities, Af e ~ 10 14 cnT 3 , redistributive collisions can interrupt 
the two-step dielectronic process itself. 

Therefore, we require dielectronic recombination from the 
metastable states as well as the ground. Secondly, we require 
final-state resolved data, i.e. we need to know the specific 
level that each two-step recombination ends-up in - the gen- 
eralized collisional-radiative population rate equations govern 
the subsequent time-evolution of these states. (Incidentally, the 
first requirement imposes an additional requirement on the sec- 
ond, namely, that we now require dielectronic recombination 
into metastable parent final-states that lie above the ioniza- 
tion limit so that the collisional-radiative modelling recovers 
the Saha-Boltzman populations and ionization fractions in the 
high-density limit.) As we have already indicated, there are 
very many possible final states accessible to the dielectronic 
recombination process which would seem to make for unac- 
ceptably large tabulations. However, it is not only impractica- 
ble, but unnecessary, to treat each final state in the same manner 
for the purposes of collisional-radiative modelling. The tech- 
niques of matrix condensation and projection reduce the effec- 
tive number of high-lying states by the progressive bundling of 



representative states and project the full influence of the high- 
lying states down onto a fully-resolved low-lying set. (It is the 
low-level set which is the focus of detailed spectral analysis.) 
With this in mind, we can tailor our dielectronic recombination 
tabulations to reflect this situation. 

Of course, many people have calculated both partial and 
total dielectronic recombination rate coefficients and it is im- 
practical to list them all here. A useful starting point is the 
compilation of total recombination rate coefficients from the 
literature by Mazzotta et al. (1998) who then used this data to 
compute the coronal ionization balance for all elements up to 
Ni. When the results of ab initio calculations are not available 
then much use, and abuse, is made of the General Formula of 
Burgess (1965). Partial dielectronic recombination rate coeffi- 
cients also abound in the literature in connection with the study 
of particular physical problems: for example, at low tempera- 
tures, where only a few autoionizing states contribute signifi- 
cantly (Nussbaumer & Storey 1983); satellite lines (e.g. Bely- 
Dubau et al. 1979); laser produced plasmas, where the electron 
density and/or charge state is high enough that the non-LTE 
populations are concentrated in a limited number of low-lying 
states (see e.g. Abdallah & Clark 1994). 

We have found it helpful for application to prepare and han- 
dle dielectronic recombination data in a hierarchy of increas- 
ing sophistication which we call baseline, level 1 and level 2. 
Baseline data are those produced using the Burgess general for- 
mula (GF), or with the techniques and state-selective programs 
associated with it. The generality of dielectronic data in use 
today are still from the GF. The background codes to the GF 
have a capability significantly beyond that of the GF, but are 
less well known. We return to these in more detail in section 
12.31 The level 1 approach was introduced in the early 'ninties' 
to support the GCR modelling of light elements such as Be, 
B and C. These species are used in iight-element strategies' 
for plasma facing wall components in fusion technology. For 
such elements, it is sufficient to use LS-coupled atomic struc- 
ture and term population modelling. The relevant metastable 
populations in this case are terms such as C 2+ (2s2p 3 P). The 
level 2 approach, which is the main purpose of the present pa- 
per series, is concerned with the need to handle medium and 
heavy species in both fusion and astrophysics and to handle 
more extreme environments. In this approach, we aim to work 
with levels rather than terms and to use an intermediate cou- 
pling scheme, based-on the use of the Breit-Pauli Hamiltonian 
(Badnell 1997). 

There are several reasons why the level 2 approach is now 
necessary: astrophysical spectral diagnostics tend to be based 
on levels rather than terms and the competition between au- 
toionization and radiation makes it difficult to partition term- 
resolved dielectronic recombination data over levels; nuclear 
spin-orbit mixing is important even for low-Z ions, e.g. carbon, 
because only a weak mixing of LS-forbidden autoionization 
rates with LS-allowed can give rise to 'forbidden' autoioniza- 
tion rates that are comparable with the dominant radiative rates 
(see Nussbaumer & Storey 1984, Badnell 1988); as Z increases 
further then LS-forbidden radiative rates start to become sig- 
nificant; finally, dielectronic recombination via fine-structure 
transitions is completely absent in LS-coupling, giving rise to a 
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large underestimate of the low-temperature dielectronic recom- 
bination rate coefficient in some iso-electronic sequences (see 
Savin et al. 1997). 

The goal of this work is to calculate multi-configuration 
intermediate coupling dielectronic recombination rate coeffi- 
cients from the (ground plus) metastable levels of an ion to 
all possible final states, resolved by level, and/or bundling, ap- 
propriate for generalized collisional-radiative modelling. We 
will cover elements applicable to astrophysics and magnetic 
fusion viz. He, Li, Be, B, C, N, O, F, Ne, Na, Mg, Al, Si, P, 
S, CI, Ar, Ca, Ti, Cr, Fe, Ni, Zn, Kr, Mo and Xe. The first 
phase of the work will be the H- through Ne-like sequences. 
Level 1 LS-coupling data for many elements of these sequences 
was calculated by Badnell (1991-92, unpublished) and incor- 
porated into the atomic database part of the Atomic Data and 
Analysis Structure (ADAS) and is routinely drawn into the 
generalized collisional-radiative part of ADAS (see Summers 
2001). So, we already have a clear pathway through to the 
complete utilization of the detailed level 2 data that we will 
produce. The second phase of the work will cover the Na-like 
through Ar-like sequences, piloted initially by further {level 1) 
LS-coupling calculations to extend the 1991-92 data. Note that 
ADAS uses year numbers for the introduction of new approxi- 
mations bringing substantive contributions to the database. The 
LS-coupled work of Badnell above has the year number '93'. 
The third phase will focus on the remaining sequences of par- 
ticular elements of interest, e.g. Fe. Even with the compactifi- 
cation of the partial dielectronic recombination data, along the 
lines already indicated, full publication in a paper journal is im- 
practical, and not especially useful. So, the entire data will be 
made available via the World Wide Web (see section |5Jl. The 
organization of the dielectronic data product follows the data 
format specifications of the ADAS Project. Dielectronic data 
are assigned to the data format adf09 (Summers 2001). 

The plan for the remainder of the paper is as follows: in 
section l2~T1 we review the generalized collisional-radiative ap- 
proach encapsulated in ADAS and which influences our ap- 
proach to handling dielectronic recombination data. In section 
12.21 we describe and justify the theoretical approach that we 
take to calculate level 2 data. In section l2~3l we review in some 
detail the essence of the Burgess approach. It will be shown 
that this remains of importance for a full exploitation of the 
new work, for example, as applied to the /-redistribution of au- 
toionizing states. Also, it is necessary to assess the progress in 
the precision of new dielectronic data in comparison with the 
baseline data. In section [5] we discuss the experimental situa- 
tion for verifying dielectronic recombination data and the role 
of external fields. Some comparisons of results of our theoret- 
ical approach with high-resolution experimental results from 
storage rings are given. In section |4] we address derived data 
and we present some illustrative comparisons of GCR effective 
coefficients obtained using baseline, level 1 and level 2 data 
from dielectronic calculations. Also, we illustrate metastable- 
resolved ionization fractions. In section[5]we give more detail 
of the organization of the database and the computational im- 
plementation of its production. We finish with a short summary. 
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2. Theory 

2.1. Generalized collisional-radiative modelling 

Consider ions X +z of element X of charge state z. We sepa- 
rate the levels of X +z into metastable levels X+ z , indexed by 
Greek indices, and excited levels, indexed by Roman indices. 
The metastable levels include the ground level. We assume that 
the excited levels Xt z are populated by excitation from all lev- 
els, p and j, of X +z , by ionization from the metastable levels of 
X+ : ~ l and recombination from the metastable levels of X* z+1 . 
The dominant population densities of these ions in the plasma 
are denoted by N z p , N z ~ l and N z+l . The excited-state population 
densities, Nf, are assumed to be in quasi-static equilibrium with 
respect to the metastable populations. Thus, 

o = X c l N o - X ( 'n x ) + X ' A ; ' + X ^ ' v:: i > 

cr 7 Li v 

where Cr. are elements of the collisional-radiative matrix de- 
fined by 

N e q r A)^+N^ Hi , j>i, (2) 

where yV e is the electron density, is the electron-impact de- 
excitation rate coefficient and A^_ >j is the spontaneous radiative 
rate, both for the j —> i transition in the ion X +z . The equivalent 
expressions for upward transitions (j < i) and for i and/or j 
replaced by a metastable index follow trivially. The diagonal 
element C z denotes the loss rate coefficient from the excited 
state i and is given by 

where 

* ; -X*v (4) 

is the total ionization rate coefficient out of i. Finally, Sf~ is the 

partial ionization rate coefficient out of metastable yu and R z * 
is the partial recombination rate coefficient out of metastable 
level v, both into level i of the ion X +z . S z ~ l includes contribu- 
tions from both direct ionization and excitation-autoionization. 
R z * includes contributions from three-body, radiative and di- 
electronic recombination. 
Solving for N z , we have 

n j = -X^^-Z^ 15 ^ 1 

<r,i fi,i 

-X<r»„ ; r,;\:- 

v,i 

cr f_i v 

-X Xv , -X'v-X'' v - (5) 

cr [IV 

where X N Z , W z and R N Z are the effective populations of /' due 
to excitation, ionization and recombination from their respec- 
tive metastables, and X,I,R, FL the corresponding coefficients. It 
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is here that the connection with spectral analysis is made. The 
total emissivity in the line j — > k is given by 



(6) 



with the populations N z given by equation The correspond- 
ing photon emissivity coefficients are defined by 



X,I,Rrnnpz _ X,I,R<rr:. at 



(7) 



respectively. Thus, the contributions to the total j —> k emis- 
sivity from excitation, ionization and recombination are given 
by " 



p 

respectively. 

The dynamic metastable populations M of X z satisfy 



(8) 



zX N z 

PJ J<r 



X N Z 



+ y is^N*- 1 + y s z - 1 x N : - i \ 

P I m 7 

cr v m / 

- Y\q z .v >\v : K .v '1 

/ j I ^op P / ; " crm "mp | 
cr v m J 



mp 



= y S : N*-Y X z ^ z N: 

/ i CD:<r—>p cr / j CD:p— ><x / 
cr cr 

p. v 

+ y arn" Z ^Vf 1 - y ffS^ 1 M 
/ i CD:v->p v / j CD:p^f/ P 

V p 

+ Z 2 ^ -X^<iv -r-\- 



(9) 



which defines the generalized collisional-radiative excita- 
tion (Xcd), ionization (ScdX recombination (Rcd) and parent 
metastable cross-coupling (Qcd) rate coefficients. (We note 
that Q z p(T = initially.) This set of equations, together with 
those for N 1 ., are sufficient to solve the low-level problem - 
those levels with principal quantum number n < n c , say. In the 
absence of dielectronic recombination, or at sufficiently high 



electron densities (e.g. > 10 18 cm -3 , say, as found in laser- 
produced plasmas) then n c can be small enough for this to be a 
complete solution since all higher levels are in collisional LTE, 
with a Boltzman population distribution. However, the pres- 
ence of dielectronic recombination in low- to medium-density 
plasmas means that n c can be prohibitively large (« 500, say). 

We now consider a projection-condensation approach that 
allows for the effect of the high-level populations (n = n > n c ) 
on the low-level populations (n < n c ). We work in the bundled- 
n picture. Here the populations are grouped according to their 
parent level and principal quantum number. We assume that 
the high-level populations are in quasi-static equilibrium with 
the low-level populations and adjacent stage metastables. Thus, 
for each parent r, the high-level populations (denoted by rn) 
satisfy 



o = y c z - n z + y c z - ^n z -, 

/ l Tnjn Til / j Tll,Tll Til 

n 77' 

+ y s-'A/r' + y^- 1 ^ 1 * 

/ J TTl,fi P / j Tfl,V v 



(10) 



which is of the same form as equation Q for the low-level 
excited-states. Thus, the low-level populations satisfy equa- 
tions of the same form as (|9j- This is a full solution for all 
levels, in the bundled-n picture. It includes direct couplings 
between the low-level populations (e.g. n — > n') and indirect 
couplings via the high-level populations (e.g. n — » n — > «')■ 
However, we already have a description of the low-level prob- 
lem in the fully-resolved picture, given by equations Q and 
We can supplement these equations with the indirect cou- 
plings of the bundled-n picture, expanded over the low-level set 
using level weighting factors, oiy. This projection corresponds 
to solving equations Q and with C, S, R and Q (includ- 
ing their appearance in equation (JSJl) replaced by C, S, 'R and Q 
where C = C+ int t, S = S + ind S , <R = R+ md R and Q = Q+ ind Q, 
where 



ncW _ 
^pcr ~ "Jp,™ 



indo z 



$1, = 



cj t „' crYc z -,(C z )~l -C z - 

77.77' 



indpz+ 1 



"R z+ , 1 = a> ivn > C z -,(C z yi -R z t\ 

„. hvn Z—i vn.vn 'vn ,vn vn,v' 



and 

indgz 



pa 



> S z ~l C J "M , R 



■z 

Tn,cr ' 



(ID 

(12) 
(13) 

(14) 



Although we now have a complete solution in terms of the 
fully-resolved low-level and bundled-n high-level picture, one 
further step is of practical significance. In order to span a wide 
range of electron densities it is necessary to treat very large 
principal quantum numbers in order to reach the collision limit 
at low densities. It is not necessary to treat each n individually. 
Rather, a set of representative n-values can be used instead. If 
Nw denotes the bundled-n populations for n = n c + 1, n c + 2, . . ., 
and N= denotes a subset of them, then the two are related via 
Nw = a)-=N=, where u-= are the interpolation coefficients. 



N. R. Badnell et al.: Dielectronic recombination data for dynamic finite-density plasmas 



5 



Substituting for Nw into equations Q and yields a con- 
densed set of equations for N= and the C, S, % and Q obtained 
from the condensed set of equations are identical in form to 
those obtained from the full set of equations. 

We note here that we have made an assumption viz. that, 
following dielectronic capture, the autoionizing state is not per- 
turbed by a further collision before it either autoionizes or ra- 
diates. This is not due to a limitation of ADAS but rather a 
choice that we have made (and defined within the adf09 specifi- 
cation) so as to make the general collisional-radiative problem 
tractable over a wide range of electron densities. Working ex- 
plicitly with autoionization and radiative rates and bound and 
non-bound states rather than partial dielectronic recombina- 
tion rate coefficients and (mostly) bound states vastly increases 
the data requirements, in general. Although our collisional- 
radiative model goes over to the correct LTE limit at high elec- 
tron densities, there is a density range (> 10 16 cm 4 , found 
in laser-produced plasmas) where the levels of spectroscopic 
interest have non-LTE populations that are influenced by non- 
LTE populations of autoionizing levels that are themselves col- 
lisionally redistributed. We describe an approximate solution in 
section l2~3l below. 

We are now in a position to spell out our requirements of 
recombination data: 

(i) We require recombination data from all metastable levels, 
not just the ground. 

(ii) We require recombination data into particular final states. 

(iii) We require the final-state to be level-resolved for n < n c 
and parent-level-resolved bundled-n for n> n c . 

(iv) Parent metastable cross-coupling means that we require re- 
combination into metastable autoionizing final states. 

(v) We need only produce data for a representative set of n. 
The ADAS adf09 data specification incorporates all of these 
requirements. Finally, we note that use of total zero-density 
ground-state recombination rate coefficients is quite unsafe for 
the collisional-radiative modelling of dynamic finite-density 
plasmas. 

2.2. Dielectronic recombination rate coefficient 
modelling 

We have already noted that the partial recombination rate coef- 
ficient (R^ 1 ) includes contributions from three-body, radiative 
and dielectronic recombination. In ADAS, three-body recom- 
bination rate coefficients are obtained from electron-impact 
ionization rate coefficients, via detailed balance. This also en- 
sures that the correct Saha-Boltzman limit is reached at high 
electron densities. Since three-body recombination is separate 
from dielectronic and radiative recombination, it is not nec- 
essary to consider if further. However, quantum mechanically, 
dielectronic and radiative recombination are indistinguishable 
processes which interfere with each other. In practice (see 
Pindzola et al. 1992), this interference is a very small effect 
and can safely be neglected for our purposes. This is the inde- 
pendent processes approximation whereby dielectronic and ra- 
diative recombination can be considered separately and is the 
approach taken by the database aspect of ADAS. Separate data 



files exist for dielectronic (adf09) and radiative recombination 
(adfO S)and they can be updated independently. Our focus is di- 
electronic recombination. Details of the ADAS data status for 
radiative recombination can be found in Summers (2001). 

In the isolated resonance approximation, the partial di- 
electronic recombination rate coefficient or: +1 from an initial 

IV 

metastable state v into a resolved final state i of an ion X +z is 
given by 

.3/2 



,c+l 



k B T e 



>v,EJ A )- 



(15) 



-j^h £jm,l* 1 j^>m,E c l 

where a)j is the statistical weight of the (A^+ l)-electron doubly- 
excited resonance state j, u> v is the statistical weight of the N- 
electron target state and the autoionization (A a ) and radiative 
(A r ) rates are in inverse seconds. Here, E c is the energy of the 
continuum electron (with angular momentum 0, which is fixed 
by the position of the resonances, and /h is the ionization po- 
tential energy of the hydrogen atom (both in the same units of 
energy), is the Boltzman constant, T e the electron tempera- 
ture and {Ancffi 12 = 6.601 1 x 10~ 24 cm 3 . The effect of interact- 
ing resonances on dielectronic recombination has been investi- 
gated by Pindzola et al. (1992) and can safely be neglected, at 
least in the absence of external electric and magnetic fields (see 
section 3 below). While autoionization rates can be determined 
(within the isolated resonance approximation) via the fitting 
of resonances calculated in a close-coupling approximation, 
or via the extrapolation of threshold close-coupling collision 
strengths using the correspondence principle, it is usual now 
to introduce a further approximation - that of using distorted 
waves, i.e. the autoionization rates are calculated via perturba- 
tion theory using the Golden Rule (Dirac 1930). This is the only 
approximation that we have made so far that may need to be re- 
considered in certain cases. In low-charge ions, a perturbative 
distorted wave calculation may give inaccurate autoionization 
rates compared to those calculated in a close-coupling approx- 
imation. However, this only has a direct effect on the partial di- 
electronic recombination rate if the autoionization rates do not 
'cancel-out' between the numerator and denominator of Jl 5i 
- typically, autoionization rates are orders of magnitude larger 
than radiative rates. 

One could obtain (some) partial dielectronic recombination 
data from an /^-matrix photoionization calculation, on making 
use of detailed balance, either in the absence of radiation damp- 
ing (Nahar & Pradhan 1994) or with its inclusion (Robicheaux 
et al. 1995, Zhang et al. 1999). (One must take care not to 
double count the radiative recombination contribution in the 
modelling now.) However, this cannot provide us with a com- 
plete set of partial recombination rate coefficients since it is 
only possible to compute photoionization from (i.e. photore- 
combination to) a relatively low-lying set of states - up to 
n ^ 10, say. Total recombination rates are obtained by sup- 
plementing the photorecombination data with high-n 'close- 
coupling' dielectronic recombination rate coefficients calcu- 
lated using Bell & Seaton (1985) or Hickman's (1984) ap- 
proach. This is based on the radiative-loss term from a unitary 
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S -matrix and does not, and cannot, resolve recombination into a 
particular final state, which is essential for collisional-radiative 
modelling. Indeed, even when summed-over all final states, er- 
rors can still result. This has been demonstrated explicitly by 
Gorczyca et al. (2002) in the case of Fe 17+ . They found that 
only the IPIRDW approach could reproduce the measured di- 
electronic recombination cross section of Savin et al. (1997, 
1999) for high Rydberg states. Thus, our initial goal is to gen- 
erate complete data sets within the independent processes and 
isolated resonance using distorted waves (IPIRDW) approxi- 
mation. Subsequently, selective upgrades from /^-matrix data 
may be made via, for example, the RmaX network which can 
be viewed as a progression of the Iron Project (Hummer et al. 
1993) and which is focusing on X-ray transitions - see, for ex- 
ample, Ballance et al. (2001). We note that while the Opacity 
Project (Seaton 1987) calculated a large amount of photoion- 
ization data, which in principle could be used for recombina- 
tion (via detailed balance), unfortunately, only total photoion- 
ization cross sections were archived, i.e. summed-over the final 
electron continuum, and so it is impossible to apply detailed 
balance and so it cannot be used as a source of recombination 
data. 

We use the code autostructure (Badnell 1986, Badnell & 
Pindzola 1989, Badnell 1997) to calculate multi-configuration 
intermediate coupling energy levels and rates within the 
IPIRDW approximation. The code can make use both of non- 
relativistic and semi-relativistic wavefunctions (Pindzola and 
Badnell 1990). The low-n problem is no different from the one 
of computing atomic structure. The high-n problem requires 
some discussion. The mean radius of a Rydberg orbital scales 
as n 2 and so it rapidly becomes impossible to calculate an ex- 
plicit bound orbital (for n > 20, say) and some approximation 
must be made. We note that (Seaton 1983) 

/ w 3\ 3 / 2 

Jim I — I P„,(r) = F kl (r)\ k=0 , (16) 

where the bound orbitals are normalized to unity and the con- 
tinuum orbitals to n6(k - k'), here k 2 = ^(Ry). The approach 
taken in autostructure is to make use of d!6i at finite n, i.e. to 
approximate the bound orbital by a suitably normalized zero- 
energy continuum orbital, for n > 15+/ 2 /4 (evaluated in integer 
arithmetic). A further refinement is to evaluate the (true) con- 
tinuum orbital for the incident electron at E c + z 2 /n 2 Rydbergs 
(instead of E c ) so as to maintain the same transition energy. 
(In the Bethe approximation, the free-free dipole acceleration 
integral is proportional to As 2 times a slowly-varying-with- 
energy dipole-length integral, where As is the transition en- 
ergy.) Actually, the true continuum orbital is calculated at about 
15 energies per I and the one- and two-body bound-free inte- 
grals are interpolated at the required energy, which is given by 
energy conservation. (The use of a zero-energy continuum or- 
bital means that long-range free-free integrals arise and these 
are treated using the techniques of distorted wave scattering 
theory, see Badnell 1983.) For each nl, autostructure reforms 
both the N- and (N+ l)-electron Hamiltonians and diagonalizes 
them (separately) to reform the rates. Only for H-like ions have 
we found it necessary to treat all / at the same time, and treat 
only each n separately. Typically, each n is calculated explicitly 



until no new continua can open-up and then only a representa- 
tive set of n, up to n — 999, is used. This approach avoids the 
extrapolation of low-n autoionization rates or, even worse, par- 
tial dielectronic recombination rate coefficients to high-n. Use 
of dl6> is, in effect, an interpolation since it is exact in the limit 
n — > oo and any error is bounded at the lowest-n by knowledge 
of the 'exact' result obtained from using P„i directly, rather than 
F Id. 

autostructure is implemented within ADAS as ADAS701 . 
It produces the raw autoionization and radiative rates. To pro- 
duce partial dielectronic recombination rate coefficients, ac- 
cording to the prescription of Sect. 2.1, requires further non- 
trivial organization of the raw data. In particular, radiative tran- 
sitions between highly-excited Rydberg states are computed 
hydrogenically and added-in during a 'post-processing' exer- 
cise with the code adasdr, which is implemented with ADAS 
as ADAS702. Also, observed energies for the core and par- 
ent levels are used at this stage to ensure accurate positioning 
of the resonances and, hence, accurate low-temperature rate 
coefficients, adasdr outputs directly the adf09 file for use by 
ADAS. Separate adf09 files are produced for different 'core- 
excitations' (n — > n'), e.g. 1 — » 2, 2 — » 2 and 2 — > 3 for Li-like 
ions. This enables selective upgrades of the adf09 database. 

2.3. Exploitation of the Burgess-Bethe approach 

We are concerned with how the precise calculations described 
above relate to other calculations and, in particular, to those 
commonly used in astrophysics. Our baseline calculation is 
based on the methodologies of Burgess, which represent what 
can be achieved without recourse to the detail of the above 
sections. The Burgess GF itself was in fact a functional 
fit to extended numerical calculations. The associated code, 
with extensions, we call the 'Burgess-Bethe general program' 
(BBGP). It will be shown in this sub-section how the BBGP 
can be used to obtain a working model for the /-redistribution 
of doubly-excited states and, hence, provide a correction to ac- 
curate, but unredistributed, dielectronic data so as to model the 
dynamic part of the plasma microfield. Also, our baseline, cal- 
culated using the BBGP, will allow an assessment of the typical 
error that is present in the general dielectronic modelling in as- 
trophysics to date. 

In the LS-coupled term picture, introduce a set P of par- 
ent terms y p 2s '' +l L p of energy E p relative to the ground parent 
term, indexed by p. Suppose that the excited parents are those 
with p m j„ < p < P. The metastable parents, which are the 
initial metastables for recombination and the final parents on 
which the recombining excited nZ-electron is built, are the sub- 
set 1 < p < p m i n = M. Let p' denote an initial parent with 
the incident electron denoted by k'l' . We wish to re-establish 
the expressions used by Burgess in his development and it is 
helpful to work in z-scaled dimensionless coordinates. Then, 
introducing z e s, the effective charge, the collision strength for 
a dipole excitation of y p 2s i> +x L p to y' p 2s p +l L' p , evaluated in the 
Bethe approximation, is given by 

Q.((S P L' p )k'l',(S p L p )klSL) = 48((z eff + l) 2 /z* ff ) 
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7 H \((2S + 1)(2L+ 1) 



,Ae pp ,)\ 2(25„ + l) 



/ v l 



x(2S p + 1)(2L; + l)fs„L^s„L„ \(F k >i> \p- 2 \ F kl )\ 



(17) 



in terms of the parent oscillator strength /, where /> = 
max(7, 1'). For dielectronic recombination, it is convenient to 
express this in terms of the Einstein A-coefficient for the parent 
transition and to analytically continue the collision strength to 
negative energies for the kl electron, that is as k = k/z^g — > i/n. 
Formally, we identify a band of free-electron energies dE/Iu = 
z 2 ff de// H , with the separation between n-shells 2z 2 ff / H /« 3 , and 
let 

\(Fve |p- 2 |^ / )| 2 d(6// H ) -» {nl\6){elInf\{F vv \p\P nl }\ 2 (18) 

Then, the (dielectronic) resonance-capture cross section is 
given by 



Q c {{S p L')k!l',{S p L p )klSL)A{EII n ) 



( (Zeff + If 



f6n 2 a 3 ^ 



I H \((2S + 1)(2L+ 1) 



.\L' P L p 1 
'{I V L 



2(2S P + 1) 
2 'AHS p L p S p L' p ) 



(Zeff + l) 4 



\(Fk'l> \P\ P nl)\ ■ d9) 



The inverse (Auger) rate coefficient is obtained by invoking de- 
tailed balance as 



A a ((S p L p )nlSL^(S p L' p )k'l') 



(Zeff + 1)' 



3 

2^? 



'Ae B 



xK^lpIP^)! 2 . 



A I (S p L p — > SpL'p) 

(Zeff + l) 4 



(20) 



For the generalized collisional-radiative modelling of light el- 
ement ions, it is convenient to use L-averaged doubly-excited 
levels, but still resolved by spin S, whereas the BBGP treat- 
ment uses S L-averaged levels. The corresponding Auger rates 
are then 



A\(S p L p )nlS -> (S p L' p )k'l') = 
l> \( A r (S p L p — > S p L' p ) 



( (Zeff + l) 2 



"eft' 



1 

2a* 



Ae, 



{(21 +1) 



(Zeff + l) 4 



\(F k . V \p\Pnl)\ 2 , 



(21) 



in both cases. The matching resonance-capture coefficients are 
obtained by detailed balance, or by summing and averaging 
over the resolved expression for Q c given above. 

Turning to the radiative decay of the doubly-excited reso- 
nant states in the LS -resolved picture: the spontaneous emis- 
sion coefficient, with a passive spectator in the nl shell, is given 
by 



A\(S p L p )nlSL -> (S p L' p )nlSL!) = (2L P + 1)(2L' + 1) 
x< /' ,V A A T (S p L p -» S p L' p ) , (22) 



L L' I 



where A T {S p L p — > S p L' p ) is the parent-core spontaneous transi- 
tion probability. The L- and LS -averaged probabilities are both 



simply equal to A r {S p L p — > S p L' p ). The BBGP method exploits 
the fact that, in the dipole case, 



A 7 



l>\(F k . V \p\Pnl) 



2a\2l + l)(z eff + l) 2 z 2 



(23) 



eff 



and efficient recurrence relations are available for the genera- 
tion of hydrogenic bound-free radial integrals for all parameter 
values. It is clear that the Bethe approximation for the partial 
collision strengths can be substantially in error for < / < 2. 
Burgess introduced correction factors for the lowest partial col- 
lision strengths, based-on a comparison with more sophisti- 
cated collision calculation results that were available at the 
time. More precisely, introduce 



(24) 



con = J]n((S P L' p )k'i',(S p L p )k\ =0 i) 
I' I 

J] £l B * ±e ((S p L' p )k'l\ (S p L p )k\ =0 l) . 



The general formula for zero-density total dielectronic re- 
combination rate coefficients used a fixed-set of cori for all par- 
ent transitions. These were based-on parent 1 s — > 2p transitions 
and it is the case that the inclusion of corrections is most sig- 
nificant for parent An > 1 transitions. To exploit the BBGP 
method beyond its use for the general formula, we must estab- 
lish the population equations of the doubly-excited levels. For 
LS -averaged levels, the number densities expressed in terms of 
their deviations, b p ,„i, from Saha-Boltzmann, and referred to 
the initial parent p', are given by 



N p ,nl=N t N^i 



7ra 2 / H 



kvT e 



3/2 



, e u p,nl ■ 

/.V 



(25) 



Then, in the BBGP baseline zero-density limit, with only res- 
onant capture from the p' parent balanced by Auger breakup 
and radiative stabilization back to the same parent, we have 



b 



p,nl 



Z,>A*(p,nl^ p'k'V) 



Y,i>A\p,nl^ p'k'l')+A r (p,nl^ p',nl)j 



(26) 



In the extended BBGP program, we can also include resonant- 
capture from initial metastables other than the ground, dipole- 
allowed collisional redistribution between adjacent doubly- 
excited Z-substates by secondary ion- and electron-impact, and 
losses by 'alternate' Auger break-up and parent radiative tran- 
sition pathways. The population equations for the /-substates of 
a doubly-excited n-shell become 



2 N^nV + X 



\1'=1±1 
p-l l+l 



l'=l±l 



P-l 



N, 



p,nl 



+ X X Km^PukI' + X A P>nl^Pi 
p, = ll'=l-l pi=l 

- {Uz<f nM ^ni + N ZeS q% +l ^N p , nM 

M l+\ P 

1p 2 ,Kl'^p,nl^P2 + j/>L pun i^, p n iN pun l . (27) 



p 2 = l l'=l-1 



p,=p+l 
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These equations may be solved progressively downwards 
through the levels built on excited parents, terminating with 
levels built on the 'final' ground and metastable parents. The 
calculations yield state-selective dielectronic recombination 
coefficients to levels built on each metastable parent, together 
with Auger rates, from levels built on metastables above the 
ground parent, to lower metastable parents. The results from 
the above solution, at zero density, we call baseline data. 
We note some details of the implementation: 

2.3.1 . Collisional rates 

Ion and electron Z-redistributive cross sections are evaluated 
following the method of Pengelly & Seaton (1964). This is to 
be viewed as a very simplified treatment of the dynamic part of 
the plasma microfield associated with ions which move closer 
than their neighbours to the target. The quasi-static part of the 
ion microfield, which is approximately 1/3 of the total, is ig- 
nored at this level of analysis (see section[3]for further discus- 
sion of field effects in an experimental context). 



are critical to the 



2.3.2. Energy levels 

The energy differences E p>n is - £p, n /±i s 
problem of doubly-excited-state redistribution. They are small, 
tending to zero at large /. Thus, the redistributive cross-sections 
are very large, remaining finite in the degenerate level limit 
only because of the finite radiative lifetime of the target (which 
is short for resonant states) or through screening of the projec- 
tile by nearest neighbours. For the / = 0, 1,2 waves of the n- 
shell spectator electron, quantum defect expansions of the form 
Hp,ni s — a o + ai/n 2 are usually available. For the higher l(> 2), 
the large number of energies required can be estimated more 
economically from the dipole polarizabilities, ap^ 1 , of the par- 
ents, following Edlen(1964), as 



E P ,„i = (zi/n) 2 [l+5.23504- 4 (zi/«) 2 (»/a+l/2) + 
3/4)] + c%\zxln)\3n z -1(1 + 1))/ 
(2n(l - 1/2)1(1 + l/2)(l + l)(l + 3/2)) . 



(28) 



Small-scale auto structure runs have been used to prepare 
these data for the baseline. We use two-term quantum defect 
expansions fitted at n — > oo and n — 10. 

2.3.3. Bethe correction factors and radiative transition 
probabilities 

For our comparisons, it is important that the inputs for the 
BBGP baseline calculations are consistent with the autostruc- 
ture calculations for the level 1 and level! results. Small-scale 
autostructure runs have been used to prepare these inputs 
for the baseline. Low partial-wave collision strengths are now 
widely available (e.g. from autostructure), such that the cor/ 
can be prepared fairly easily. Within the ADAS Project, effi- 
cient subroutines for bound-free integrals and complete BBGP 
calculations are available together with sets of cor/ for the prin- 
cipal types of transitions. They are available to those for whom 
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p. = 1 S; p„ = 1 S : doublet 

p„ = 1 S; p = 3 P : doublet 
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p = 3 P; p„ = 3 P : quartet 

Pp = 3 P; p„ = 1 S : doublet 



T| 1 

Quenched • 




Fig. 1. The graphs contrast a^ ot \p p — » pi — > p a ) vs T e for the 
dielectronic recombination of 4+ ions in various approxima- 
tions at zero density, (a) p p = 4+ (2s 2 1 S), /v= 3+ (2s 2 2p 2 P): 
An — and An = 1 intermediate parents included. The 
Burgess GF includes only the dipole 2s 2 'S - 2s2p 'P and 
2s 2 'S - 2s3p 'P transitions. BBGP incorporates specific cori 
correction factors, energies which differ from those implicit in 
the GF, and the 3p - 3s alternative Auger channel. The level 
1 and level 2 results include all allowed and non-allowed par- 
ent transitions within the n — 2 and n — 3 complexes. Note 
the low temperature extension which cannot be modelled with 
the GF and BBGP. The correct distinction and positioning of 
the key lowest resonances are possible only at level 2. (b) 
p p = 4+ (2s 2 >S) and 4+ (2s2p 3 P), /v=0 3+ (2s 2 2p 2 P) and 
3+ (2s2p 2 4 P): level 1 results separated by spin-system and fi- 
nal parent. Recombination from-and-to metastables cannot be 
handled by the GF. Excited-states built on the 2s2p 3 P par- 
ent have a spin-change autoionization pathway. The level 1 
metastable-resolved totals do not include this loss. Within an 
LS-coupled GCR picture using level 1 data, spin-breakdown 
Auger data is included explicitly in the GCR calculations for 
the correct linking of systems built on the 2s 2 'S and 2s2p 3 P 
parents. The relevant final-parent-changing Auger data is in- 
cluded in the adf09 data file specification. For comparison with 
simpler treatments, totals including quenching of n-shells > 4 

Knilt rvn tti£» 7q7ii n^rpnt art* ^lcr» cVir»\n/n 
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Fig. 2. 4+ recombination, (a) Partial n-shell recombination co- 
efficients: initial recombining metastable 2s 2 'S, intermediate 
excited parent p,, final states (2s 2 'S)n with /?, = 2s2p 'P and 
2s3p 'P. r e = 1.6 x 10 6 K, N e = N p = 0. BBGP (gf) indicates 
use of standard cor/ matching the GF, BBGP (exact) indicates use 
of specific cor,. For the An = case, BBGP (gf) and BBGP (exact) 
are superposed. An = 1 cases: 1. no alternative Auger chan- 
nels, GF corf, 2. no alternative Auger channels, specific cor,; 

3. 2s3p 'P - 2s3s 'S alternative Auger channel, specific corf, 

4. as before and 2s3p 'P - 2s2p 'P alternative Auger channel, 
specific cori. (b) b PiJl \ factors for doubly-excited states of 3+ 
relative to 4+ 2s 2 'S for p, = 2s2p 'P and 2s3p 1 P, n = 20, 
T e = 10 6 K, N e = N p , Zeff = 1 . Cases: 1. N e = 10 10 cm" 3 ; 2. 
N e = 10 12 cm" 3 ; 3. N e = 10 13 cm" 3 ; 4. N e = 10 14 cm" 3 ; 5. 
A^ e = 10 cm -3 . Note the alternative Auger channel reduction 
for the pi = 2s3p 'P graphs. 



direct utilization of level 1 and level 2 data in full GCR mod- 
elling is not an option. 

The results presented in this overview paper are illustrative 
only. Fig. ^ contrasts zero-density total dielectronic recombi- 
nation coefficients (a 1 ? 1 ) calculated in the GF and the BBGP 
baseline approximations with those of the level 1 and level 2 
computations reported here. In particular, we note that level 1 
and 2 data are required to describe the recombination at low 
temperatures and that the level 2 data provides a noticeable re- 
finement over the level 1 results. Figs.|2t,b illustrate the partial 
recombination into n-shells and the population structure of the 



/-subshells of a representative doubly-excited n-shell. Fig. |5Ji 
shows the very good convergence of BBGP to level 1 data with 
increasing completeness of alternate Auger pathways. Fig.fJJ) 
shows the effects of collisional redistribution at finite-density. 
A ratio of the sum over /-substates at a given density to that 
at zero density yields a BBGP finite-density adjustment factor 
of the total n-shell capture at zero density. The consistency be- 
tween the BBGP, level 1, and level 2 approaches allows us to 
use this adjustment factor on the level 1 and level 2 data. In ad- 
vanced generalized collisional-radiative modelling, the BBGP 
finite-density redistributive code acts as an interface between 
the extraction of state-selective zero-density dielectronic data 
from the ADAS adf09 database and its entry into the GCR pop- 
ulation codes, corrected for finite-density doubly-excited state 
redistribution. Note also that routine semi-automatic compar- 
isons as, illustrated here, provide the theoretical uncertainty es- 
timate with which we can tag each dielectronic datum. 

3. Experimental validation of dielectronic 
recombination data and the role of fields 

The last decade has seen an enormous amount of experimen- 
tal activity in the area of dielectronic recombination. In par- 
ticular, heavy-ion storage rings coupled with electron-coolers 
have provided a wealth of data for partial dielectronic recom- 
bination. (The partial here is by the intermediate resonance 
state rather than the final state.) Most data is of the 'bundled- 
n' form, but some /-resolution is possible for very low-lying 
states. The iso-electronic sequences studied range almost ex- 
clusively from H-like through to Na-like and nuclear charges 
have ranged between Z=2 and 92. In all cases, in the absence 
of external fields, there is rarely any significant disagreement 
with theory, i.e. outside of the experimental uncertainty. A few 
of the more recent, typical, comparisons between experiment 
and the results of IPIRDW calculations include: Bohm et al. 
(2002), Savin et al. (2002a,b, ), Brandau et al. (2002). In Fig-EJ 
we show representative comparisons of dielectronic recombi- 
nation data for O s+ + e~ — >0 4+ and Fe 18+ + e~ — >Fe 17+ , cal- 
culated in the IPIRDW approximation with autostructure and 
which illustrates the level of accuracy that can be expected of 
the theoretical data. 

One major area of uncertainty is the role of external fields 
on dielectronic recombination, and it is this more than any- 
thing that renders pointless efforts to compute (zero-density) 
field-free data to an accuracy of better than ^20% , say. It has 
long been known that the high Rydberg states that frequently 
dominate the dielectronic recombination process can be Stark- 
mixed by weak electric fields (Burgess & Summers 1969), in 
particular the plasma microfield (Jacobs et al. 1976), and so 
increase the partial rate coefficients by factors of 2, or 3, or 
more, over a wide range of n. Recently, the picture has been 
further complicated by the discovery that magnetic fields, when 
crossed with an electric field, strongly affect the electric field 
enhancement - by reducing it in most cases (see Robicheaux et 
al. 1997, Bartsch et al. 1999, Schippers et al. 2000 and Bohm 
et al. 2001). While this suppression of the electric field en- 
hancement is advantageous towards the use of field-free di- 
electronic recombination data, it is disadvantageous in terms of 
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Fig. 3. (a) A theoretical reconstruction is shown of the ob- 
served dielectronic recombination resonances for the light ion 



O 



5+ 



+ e 



>0 4+ . The measurement is from the CRYRING 



heavy-ion storage ring (Bohm et al. 2002). Note that the di- 
electronic recombination data from the methods described in 
section I2.2I have been convoluted with the experimental ve- 
locity distribution, which has the units of a rate coefficient. 
The three pairs of curves correspond to three different beam 
energies which results-in three different cut-offs in the max- 
imum principal quantum number detected. The experimental 
and theoretical data agree quite nicely, (b) The Maxwellian 
dielectronic recombination rate coefficient for Fe I8+ + e~ — > 
Fe 17+ is shown. The measurement (solid curve) is from the 
Heidelberg storage ring (Savin et al. 2002a) with an experi- 
mental error assessed as < 20%. The theoretical IPIRDW re- 
sults include the 1 — > 2, 2 — > 2 and 2 — * 3 core excitations 
in the multi-configuration Breit-Pauli (MCBP, dotted curve) 
and Dirac-Fock data (MCDF, dashed curve) with n max = oo. 
The worst deviation from experiments is < 30% (for the 3/3Z' 
resonances) with a typical uncertainty < 20% for the direct 
state-selective coefficients to individual levels and for the total 
dielectronic rate coefficient. There is excellent agreement be- 
tween the Breit-Pauli and Dirac-Fock results. The Breit-Pauli 
approach is used for our mass data production. 

trying to compute field-dependent data for plasma modelling. 
Previously, it appeared that a reasonable approach would be 
to use the values of the plasma microfield (which in turn de- 
pends on the plasma density) for the electric field strength for 



use in the generation of field dependent (i.e. density depen- 
dent) data as input to collisional-radiative modelling. This in it- 
self ignored any further (e.g. external) electric fields that might 
be present in the plasma environment, beyond the plasma mi- 
crofield. The recognition of the importance of magnetic fields 
as well makes a comprehensive solution to dielectronic recom- 
bination in a plasma a distant goal and partial data accurate 
to ^20% as meaningful as necessary. Furthermore, field en- 
hancement is sensitive to interacting resonances as well (see 
Robicheaux et al. 1998) unlike the field-free case. We do note 
again that high Rydberg states in a finite density plasma are 
brought into LTE by (electron) collisions. A preliminary study 
by Badnell et al. (1993) showed that the effect of the plasma mi- 
crofield on the density-dependent effective recombination rate 
coefficient was suppressed by collisions driving high Rydberg 
states into LTE - larger values for the microfield, which lead to 
larger enhancements of the zero-density rate coefficient, corre- 
sponds to denser plasmas for which collisions drive more states 
into LTE. 

4. Comparisons and assessments of the derived 
data 



As described in section 12.11 the emissivity and generalized 
collisional-radiative (GCR) coefficients depend inter alia on 
the fundamental dielectronic cross section data. There are two 
issues of concern in assessing these derived theoretical data, 
viz. the relative contributions to the effective coefficients com- 
ing from many different direct and indirect pathways and, sec- 
ondly, estimation of the uncertainty in the theoretical data, 
which may be treated as a 'working error' in the interpretation 
of spectral observations from plasmas. 

For the effective photon emissivity coefficients (CP£Cs), it 
is firstly to be noted that the relative importance of the contri- 
bution from excitation ( X 3 3 £C) and recombination ( R T£C) is 
directly proportional to the ionization balance fractional abun- 
dances of the (metastable) 'driver' populations. The recombi- 
nation part is most significant in transiently recombining plas- 
mas and it is on this part only that we focus here. The parti- 
tioning of the collisional-radiative matrix described in section 
12. H allows us to contrast the direct capture, capture coming via 
the complete set of resolved low-levels and capture via the high 
bundled-n quantum shells, which are treated by projection. The 
relative contributions depend differentially on density since the 
projection part is suppressed selectively at higher densities. 
Also, electron temperature and the recombining ion charge in- 
fluence the relative importance of the dielectronic and radiative 
recombination contributions and the role of the more highly- 
excited levels. In Fig. 0] we show the main effects with some 
illustrative results from ADAS for the C in 2s2p 3 P - 2s3d 3 D 
multiplet at 459.6 A and the Ne vn 2s2p 3 P - 2s3d 3 D multiplet 
at 106.1 A. 

Fig. [5] illustrates the main features of the generalized 
collisional-radiative recombination coefficients for 4+ + 
e~ — >0 3+ . In the GCR term-resolved picture for light elements, 
there are four coefficients associated with the pairings of the 
2s 2 'S & 2s2p 'S and 2s 2 2p 2 P & 2s2p 2 4 P metastable terms 
in the recombining and recombined systems, respectively. As 
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Fig. 4. (a) Temperature and density dependence of the GCR 
R ?£C Z for the C m 2s2p 3 P - 2s3d 3 D multiplet at 459.6 A 
driven by the C 3+ (2s 2 S) recombining ion. (b) R J'£C Z for the 
C m 2s2p 3 P - 2s3d 3 D multiplet at 459.6 A and Ne vn 
2s2p 3 P-2s3d 3 D multiplet at 106.1 A, respectively. Additional 
curves contrast the corresponding CPSC^s driven by the ground 
states of the recombined ions at p = 10 12 cm 4 . For compari- 
son between isoelectronic systems, it is convenient to use the 
scaled electron temperature 6 = T e /z\ and scaled electron den- 
sity p = N e /z[, where z\ is the recombining ion charge (=3 and 
7, respectively). Since the upper level is in the excited n — 3 
shell, the cascading (projected) influence of higher levels is 
larger than for spectrum lines originating in the n — 2 shell. 
The low temperature behaviour is that of the radiative recom- 
bination process at low density, but rises to that of the collective 
three-body process at high density. 



the radiative and three-body processes are included, the low 




Fig. 5. (a) 4+ + e~ — >0 3+ term-resolved GCR recombination 
coefficients that are driven from the 2s 2 : S ground metastable 
of the 4+ ion as a function of electron temperature for a num- 
ber of electron densities. The scaled electron temperature, 6, 
and scaled electron density, p, are used, (b) 4+ + e~ — >0 3+ 
term-resolved GCR recombination coefficient driven from the 
2s2p 3 P metastable of 4+ . The suppression of the coefficients 
compared with those of (a) is due to spin-breakdown alternative 
Auger channels. These must be included even in term-resolved 
GCR modelling. 



temperature and high density behaviours reflect these contribu- 
tions. The finite-density suppression of the coefficient for the 
ground parent case and the effect of alterative Auger channels 
are both pronounced and also depend on the ion charge. These 
effects require GCR modelling. The simpler stage-to-stage pic- 
ture introduces a significant and, generally, unquantifiable er- 
ror. It is to be noted that the intermediate-coupled dielectronic 
recombination data of this project also sustains production of 
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fine-structure-resolved metastable GCR coefficients appropri- 
ate to medium- and heavy-weight elements. 

The present paper's main concern is with ensuring the qual- 
ity and completeness of dielectronic data for plasma modelling 
and not with all the consequential modelling of populations and 
ion distributions in plasmas. There are, however, two points 
to draw attention to. Firstly, it is well known that dielectronic 
recombination shifts equilibrium ionization balance fraction 
curves to higher temperatures. This is most pronounced for the 
ions with one or two electrons outside of closed shells and pro- 
duces a characteristic 'piling-up' of these stages. It is also these 
ion fractions which show most markedly the effect of finite 
density reduction of the collisional-radiative coefficients. This 
cannot be ignored for moderately ionized systems in plasmas 
with N e >10 10 ctrT 3 . These effects are shown in Fig.[6]for oxy- 
gen. Secondly, most plasma transport models work only with 
whole ionization stage populations. The present work, how- 
ever, sustains the metastable-resolved GCR picture. GCR co- 
efficients and ionization balance fractional abundances must be 
bundled back to the ionization stage for such models, at the 
expense of precision. Fig. [6j> illustrates the resolved picture 
for the beryllium-like ionization stage of oxygen. The sim- 
plest bundling strategy imposes equilibrium fractions on the 
metastable populations relative to the ground, as is used for the 
stage-to-stage fractional abundances in Fig.|(J- 

Modern good practice requires an estimate of uncertainty 
in derived theoretical data so that meaningful deductions can 
be drawn from the comparison with observations. It is unfor- 
tunately the case that most theoretical dielectronic data has no 
error associated with it. Because of the relative complexity of 
dielectronic recombination and the many contributions, agree- 
ments between different theories and with observations some- 
times appear fortuitously and do not reflect the underlying re- 
liability (see Savin et al. 2002a). For the present derived GCR 
coefficients and R CP£C, we outline our approach to procuring a 
relevant 'working error'. 

In the ADAS project (Summers et al. 2002), a distinction 
is made between 'locked' parameters, as distinct from 'search' 
parameters, in the optimized fitting of models to observations. 
Search parameters return a fit uncertainty or confidence level, 
the locked parameters must carry an error with them. An effec- 
tive rate coefficient is such a locked parameter. Its uncertainty, 
called the cumulative statistical error, is computed from the er- 
rors of the fundamental reaction rates as follows: Monte Carlo 
samples are made of all the individual reaction coefficients, 
within their (assumed) independent Gaussian uncertainty dis- 
tributions, and the derived coefficient calculated. The process 
is repeated many times until statistics are built up. The accu- 
mulated results are fitted with a Gaussian variance. 

The key issue then is the starting point of uncertainties 
in the fundamental component dielectronic coefficients. The 
BBGP codes described in section l2~3l have been arranged to 
generate adf09 baseline files. Such a file is differenced with 
the matching level 1 file to provide an error estimate for the 
baseline values and may be stored in a .err file exactly parallel- 
ing the naming of the actual .dat file. In like manner, the level 
1 file may be differenced with the level 2 file (averaged-over 
fine-structure) to provide the .err file for level 1. We treat this 



Oxygen coronal equilibrium 



GCR (zero density) 
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4+ coronal equilibrium 



1 4 CITT 
1 8 CITT 
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Fig. 6. (a) Behaviour of the ionization balance fractional abun- 
dances for oxygen as a function of electron temperature and 
density. The balance is also shown using the GF for the dielec- 
tronic contributions. The same effective ionization rate coeffi- 
cients were used in all three curve sets, but are not considered 
further here. The (level 1) GCR results are shown as a stage- 
to-stage balance, but originating from a true GCR metastable- 
resolved calculation. The metastable fractions are combined by 
weighting with their equilibrium fractions as determined by a 
low-level population balance. Note the potential confusion be- 
tween differences due to the use of a low precision zero-density 
dielectronic calculation, such as the GF, and those due to finite- 
density effects. A more complete and sophisticated approach to 
such metastables (extended also to ionization stages) is called 
'flexible partitioning' and will be the subject of a separate 
work, (b) Beryllium-like ionization stage, 4+ , fractional abun- 
dances in the metastable term-resolved GCR picture. Curves 
are shown, therefore, for both the 2s 2 'S and 2s2p 3 P ground 
and metastable terms. (The metastable curves are completely 
supressed at the lowest two densities.) 
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3+ 2s2p 2 2 P-* 2s 2 2p 2 P 




Fig. 7. Upper and lower cumulative statistical (+1 standard de- 
viation) surfaces for the R ?£e z for the O iv 2s2p 2 2 P-2s 2 2p 2 P 
multiplet at 554.4 A. Note that projection has a substantial in- 
fluence on the line emissivity and so its influence on the cumu- 
lative propagated error is significant. Projection and its error 
contribution is not included in the graph shown. 78.C Z data are 
archived in ADAS data format adfl5. The propagated 'locked 
parameter' errors are available in .err files paralleling the .dat 
files for use in plasma modelling. 



also as a conservative error for level 2. It is emphasized that 
this is not a confident absolute error, but a (hopefully) helpful 
appraisal of the theoretical data. It is most appropriate for the 
n- and nl- shell bundled data. The experimental comparisons 
of the type discussed in section [^indicate that a minimum un- 
certainty ^20% is appropriate for the term and level selective 
dielectronic data. Use of such .err files is not yet a common 
practice and its handling within a projection matrix framework 
is complex. R 3 £C < error surfaces are shown in Fig. fusing the 
ADAS procedure, but propagating error only from the state- 
selective part. The full handling of error will be treated in a 
separate paper. 

5. Structure and access to the database 

The complete set of dielectronic recombination data (both 'par- 
tial' and 'total') will be publically available as adf09 files from 
the Controlled Fusion Atomic Data Center at the Oak Ridge 
National Laboratory, USA (http://www-cfadc.phys.ornl.gov/!. 
These data files are simple ascii text in a formatted organiza- 
tion. The layout differs slightly between the level 1 LS-coupled 
and level 2 intermediate-coupled forms. Fig. [8] summarizes 
the intermediate-coupled form. The document appxa_09.pdf 
of Appendix A of the ADAS user manual (available at: 
http://adas.phys.strath.ac.uki provides the detailed description. 
This includes a summary of the sublibraries and their current 
status, content of the data lines and the meanings of all param- 
eters, together with some samples of the format. An example of 



Specify 

element ion and coupling 



Recombining ion: 

level indexing and energies 



Recombined ion: 

resolved level indexing and energies 



Recombined ion: 

nl-shell indexing and Auger rates 



Recombined ion: 

n-shell indexing and Auger rates 



Cycle over initial parent levels 



Resolved level DR coeffts. Vs Te 



Cycle over final parent levels 



nl-shell DR coeffts. Vs Te 



n-shell DR coeffts. Vs Te 



Fig. 8. Organization of data within the adf09 format. adf09 
specifies both LS- and intermediate-coupled data organiza- 
tions. In the LS-coupled case, the coefficients span resolved 
terms with valence electron up to n — 7; « 40 representative 
n-shells up to n = 999. In the intermediate-coupled case, the 
coefficients span resolved levels with valence electron up to 
n = 7; all nZ-shells up to n = 10; ~ 40 representative n-shells 
up to n — 999. The coefficients are tabulated at 19 scaled tem- 
peratures spanning from 10 - 10 7 K. Auger rates for ioniza- 
tion to alternate metastable parents for the set of nZ-shell and n- 
shell spectators built on each parent-metastable are included for 
model completeness. The detailed specification is in Appendix 
A of the ADAS User's Manual (Summers 2001). 



a pathway to a member is /../adf09/jc00#li/jc00#lijie7ic23.dat 
which distinguishes the producer initials ' jc\ year number 
'00', recombining ion iso-electronic sequence 'IV, element 
'ne', coupling He' and parent n — 2 — > n — 3 transition group 
'23'. 

We have found it convenient to archive also the driver 
data sets, which initiate the level 1 and level 2 calculations, as 
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Drivers 
for Autostr. 



Drivers 
for Autostr. 
postproc. 

ADF28 



ADAS701 



Auto- 
structure 




yy 

ADAS702 

Post-process 
for 

dielectronic 
rate coeffts. 



ADAS703 

Post-process 
for doubly- 
excited 
populations 



ADAS410 

Extract and 

display 
dielectronic 
data 






ADAS212 

Map state 
selective 
diel. coeffts. 
to specific 
ion files 


►j 







Specific 
ion file 
for sat. 
line calc. 

ADF04 



Specific 
ion file 
forGCR 
calc. 

ADF04 



Fig. 9. ADAS code organization for production of the com- 
plete set of dielectronic data. ADAS701 is the autostructure 
code. The adasdr post-processor code, ADAS702, prepares 
the level-resolved, bundled-n/ and bundled-n partial recombi- 
nation coefficient data according to the adf09 format specifi- 
cation (Summers 2001). The code ADAS212 maps the level- 
resolved data onto the specific ion file class adf04. ADAS703 is 
an additional post-processor for dielectronic satellite line mod- 
elling and will be the subject of a separate paper. In a sepa- 
rate code chain (not shown here), ADAS807 prepares cross- 
referencing files to the bundled-n/ and bundled-n data which 
are required for the very high-level population calculations and 
the evaluation of the projection matrices by ADAS204. The 
fully-configured adf04 files, together with the projection matri- 
ces, are processed by ADAS208 which delivers the final gen- 
eralized collisional-radiative (GCR) coefficients and effective 
emission coefficients. 



other ADAS data formats (adf27 and adf28). These have path- 
ways which parallel adf09. Various ADAS codes execute the 
primary calculations and the subsequent collisional-radiative 
modelling. The flow of calculation is summarized in Fig. [5] 



6. Summary 

We have described the goals and methodology behind a pro- 
gramme to calculate a comprehensive database of dielectronic 
recombination data for the collisional-radiative modelling of 
dynamic finite-density plasmas and illustrated its use in such 
environments. The first phase of the program covering H- 
through Ne-like sequences is under way and illustrative results, 
comparisons, and total (zero-density) rates will be the subject 



of a series of papers to be submitted to A&A in the near future, 
e.g. O-like ions, Zatsarinny et al. (2003). 
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